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Abstract 

Intermediate statistics interpolating from Bose statistics to Fermi statistics are 
allowed in two dimensions. This is due to the particular topology of the two dimen- 
sional configuration space of identical particles, leading to non trivial braiding of 
particles around each other. One arrives at quantum many-body states with a multi- 
valued phase factor, which encodes the anyonic nature of particle windings. Bosons 
and fermions appear as two limiting cases. Gauging away the phase leads to the 
so-called anyon model, where the charge of each particle interacts "a la Aharonov- 
Bohm" with the fluxes carried by the other particles. The multivaluedness of the 
wave function has been traded off for topological interactions between ordinary 
particles. An alternative Lagrangian formulation uses a topological Chern-Simons 
term in 2+1 dimensions. Taking into account the short distance repulsion between 
particles leads to an Hamiltonian well defined in perturbation theory, where all per- 
turbative divergences have disappeared. Together with numerical and semi-classical 
studies, perturbation theory is a basic analytical tool at disposal to study the model, 
since finding the exact iV-body spectrum seems out of reach (except in the 2-body 
case which is solvable, or for particular classes of iV-body eigenstates which gener- 
alize some 2-body eigenstates). However, a simplification arises when the anyons 
are coupled to an external homogeneous magnetic field. In the case of a strong 
field, by projecting the system on its lowest Landau level (LLL, thus the LLL-anyon 
model), the anyon model becomes solvable, i.e. the classes of exact eigenstates al- 
luded to above provide for a complete interpolation from the LLL-Bose spectrum to 
the LLL-Fermi spectrum. Being a solvable model allows for an explicit knowledge 
of the equation of state and of the mean occupation number in the LLL, which do 
interpolate from the Bose to the Fermi cases. It also provides for a combinatorial 
interpretation of LLL-anyon braiding statistics in terms of occupation of single par- 
ticle states. The LLL-anyon model might also be relevant experimentally: a gas of 
electrons in a strong magnetic field is known to exhibit a quantized Hall conduc- 
tance, leading to the integer and fractional quantum Hall effects. Haldane/exclusion 
statistics, introduced to describe FQHE edge excitations, is a priori different from 



anyon statistics, since it is not denned by braiding considerations, but rather by 
counting arguments in the space of available states. However, it has been shown to 
lead to the same kind of thermodynamics as the LLL-anyon thermodynamics (or, in 
other words, the LLL-anyon model is a microscopic quantum mechanical realization 
of Haldane's statistics). The one dimensional Calogero model is also shown to have 
the same kind of thermodynamics as the LLL-anyons thermodynamics. This is not 
a coincidence: the LLL-anyon model and the Calogero model are intimately related, 
the latter being a particular limit of the former. Finally, on the purely combinatorial 
side, the minimal difference partition problem -partition of integers with minimal 
difference constraints on their parts- can also be mapped on an abstract exclusion 
statistics model with a constant one-body density of states, which is neither the 
LLL-anyon model nor the Calogero model. 
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1 Introduction 

Quantum statistics, which is concerned with quantum many-body wavefunctions of iden- 
tical particles, has a long history going back to Bose and Fermi. The concept of statistics 
originates at the classical level in the Gibbs paradox, which is solved by means of the in- 
discernability postulate for identical particles. At the quantum level, the usual reasoning 
shows that only two types of statistics can exist, bosonic or fermionic. Indeed, since 

• interchanging the positions of two identical particles can only amount to multiplying 
their 2-body wavefunction by a phase factor, 

• a double exchange puts back the particles at their original position, 

• and one usually insists on the univaluedness of the wavefunction, 

this phase factor can be only 1 (boson) or -1 (fermion). 

However, non trivial phase factor should be possible, since wavefunctions are anyway 
defined up to a phase. The configuration space of two, or more generally, N identical 
particles has to be defined cautiously PQ: denoting by C the configuration space of a single 
particle (C = R 2 for particles in the two-dimensional plane, d = 2), the configuration 
space of N particles should be of the type C n /Sn, where C x C x ... x C = C N and 
Sn is the permutation group for N identical particles. Quotienting by Sn takes into 
account the identity of the particles which implies that one cannot distinguish between 
two configurations related by an operation of the permutation group. One should also 
subtract from C N the diagonal of the configuration space Dn, i.e. any configurations 
where two or more particles coincide. The reason is, having in mind Fermi statistics, that 
the Pauli exclusion principle should be enforced in some way. A more precise argument is 
to have a valid classification of paths in the iV-particle configuration space, which would 
be ambiguous if two or more particles coinciding at some time is allowed (since they are 
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identical, did they cross each other, or did they scatter off each other ?). It follows that 
the configuration space of N identical particles should be 



c« = (i) 

Note that on this configuration space, a fermionic wavefunction is multivalued (two 
values 1 and -1), so there is no reason not to allow more general multivaluedness. Here 
come some topological arguments, which allow to distinguish between d = 2 and d > 2, 
and, as we will see later, which can be related to spin considerations. In 2 dimensions, C N 
is multiply connected and its topology is non trivial: it is not possible to shrink a path of 
a particle encircling another particle, due to the topological obstruction materialized by 
the latter. It follows that Cn is multiply connected. This is not the case in dimensions 
higher than 2, where C N is simply connected, meaning that all paths made by a particle 
can be continuously deformed into each other, i.e. one cannot distinguish the interior 
from the exterior of a closed path of a particle around other particles. 

These arguments imply that the equivalent classes of paths (first homotopy group) in 
Cn are, when d = 2, in one-to-one correspondence with the elements of the braid group 

U^Cn) = B N (2) 

whereas, when d > 2, they are in one-to-one correspondence with the elements of the 
permutation group 

MCn) = S N (3) 

The braid group generators <7j interchange the position of particle i with particle i + 
This operation can be made in an anti-clockwise manner (<jj) or a clockwise manner (er^ 1 ). 
Each braiding of N particles consists of a sequence of interchanges of pairs of neighboring 
particles via the oVs and the cr^ s, with i — 1,2, ...,N — 1. The braid group relations 
list the equivalent braiding, i.e. braiding that can be continuously deformed one into the 
other without encountering a topological obstruction 

cricr i+1 ai = a i+ i(Ti(Ti + i; o,o, o ,o, when \i — j\ > 2 (4) 

Saying that d = 2 is different from d > 2 is nothing but recognizing that Cj ^ aT 1 
when d = 2, whereas a, = o^ x when d > 2 (<Xj can be continuously deformed into &7 1 
when particles i and i+1 are not stuck in a plane). It follows that when d > 2, the braid 
generators Cj's defined by Q with the additional constraint o~j = a^ 1 are the permutation 
group generators. 

Note also that the d = 2 paradigm, <7j ^ err ; hints at an orientation of the plane, a 
hallmark of the presence of some sort of magnetic field. This point will become apparent 
in the Aharonov-Bohm formulation of the anyon model. 

The fact that C N is multiply connected when d = 2 and not when d > 2 can also 
be related to the rotation group 0(d), and thus to some spin-statistics considerations 
[2]. When d > 2, the rotation group is doubly connected, [IIi(0(d)) = Z 2 ], its universal 
covering, for example when d = 3, is SU(2), which allows for either integer or half 
integer angular momentum states, that is to say either single valued or double valued 
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Figure 1: The braid group generators and their defining relations. 

representations of the rotation group. On the other hand, when d = 2, the rotation 
group is abelian and infinitely connected [IIi(0(2)) = Z], its universal covering is the 
real line, that is to say arbitrary angular momenta are possible, and therefore multivalued 
representations. One can see here a hint about the spin-statistics connection, where 
statistics and spin are trivial (Bose-Fermi statistics, integer-half integer spin) when d > 2, 
and not when d = 2. 

Let us consider the simple one-dimensional irreducible representation of the braid 
group, which amounts to a common phase factor exp(— irca) for each generator a (and 
thus expect) to cr _1 ). It means that a non trivial phase has been associated with the 
winding of particle i around particle i + l. Higher dimensional representations (quantum 
vector states) are possible -one speaks of non abelian anyons, in that case not only a 
non trivial phase materializes during a winding, but also the direction of the vector state 
in the Hilbert space is affected- but they will not be discussed here (even though they 
might play a role in the discussion of certain FQHE fractions [3] , and, in a quite different 
perspective, in the definition of topologically protected fault-tolerant quantum computers 

HI)- 

Clearly, when d > 2, a, = a i 1 implies a = or a = 1, i.e. Bose or Fermi statistics 
(an interchange leaves the wavefunction unchanged or affected by a minus sign). 

From now on let us concentrate on d = 2 and denote the free many-body wavefunction 
of iV identical particles by i^'(fi, r%, fjv). Indeed, statistics should be defined for free 
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particles with Hamiltonian 

N 
i=l 

and special boundary conditions on the wavefunction, as in the Bose case (symmetric 
boundary condition) and the Fermi case (antisymmetric boundary condition). As already 
said, ip'(fi, r*2, fjv) is affected by a phase exp(— mot) when particles i and i + 1 are 
interchanged: one can encode this non trivial exchange property by defining 

ip'(n,r 2 , ...,r N ) = expj-iay^ 9 i j)ip(f 1 , f 2 , f N ) (6) 

i<j 

where ip(fi, f 2 , fv) is a regular wavefunction, say bosonic by convention, and % is 
the angle between the vector fj — f = fj and a fixed direction in the plane. Indeed 
interchanging i with j amounts to 9^ — > 0^ ± 7T, which altogether with the bosonic 
symmetry of ^(^1,^2, ■ ••,f/v), leads to 

ip'(ri,r 2 , ...r}, ...fj, ...,fv) = exp(=p7ra)?//(fi, f>, ...fj, ...f,-, ...,fv) (7) 

By the above bosonic convention for i/}{fi,f 2 , ...,fv), the statistical parameter a even 
(odd) integer corresponds to Bose (Fermi) statistics. It is defined modulo 2, since two 
quanta of flux can always be gauged away by a regular gauge transformation while pre- 
serving the symmetry of the wavefunctions in the Bose or Fermi systems. Indeed, (|6| can 
be interpreted as a gauge transformation. Let us compute the resulting Hamiltonian H N 
acting on ip(fi, f 2 , fv) 

N 1 

^ = E^-^)) 2 ( 8 ) 

1=1 

where 

— * 

A(f) = *d % (£ On) = aJ2 k ~^P l (9) 
k<i y 

is the statistical potential vector associated with the multivalued phase (the gauge pa- 
rameter). The free multivalued wavefunction has been traded off for a regular bosonic 
wavefunction with topological singular magnetic interactions. The statistical potential 
vector (J9]) can be viewed as the Aharonov-Bohm (A-B) potential vector that particle i 
carrying a charge e would feel due to the flux tube carried by the other particles, with 
e and related to the statistical parameter a by a = e<fi/(2it) — 4>/4>q (0o — 27r/e is the 
flux quantum in units h = 1). The resulting composite charge-flux picture is known under 
the name of anyon model j5] since it describes particles with "any" (any-on) statistics. 
Computing the field strength one obtains 

a z k A fj,- 27ra r,-, ^ ; > 

e^ A E^ = — E 5 ^') ( 10 ) 

meaning that each particle carries an infinite singular magnetic field with flux = 2%a/e. 
The gauge transformation is singular since it does not preserve the field strength (which 
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Figure 2: The two equivalent formulations of anyon statistics in terms, on the left, of a 
punctured plane and, on the right, of usual bosonic particles interacting via topological A- 
B interactions. The loop of particle i around particle j cannot be continuously deformed 
to nothing due to the topological obstruction materialized by the puncture at the location 
of particle j. 



vanishes in the multivalued gauge and is singular in the regular gauge). This is due to 
the singular behavior of the gauge parameter 6^ when particle i come close to particle j, 
thus the singular Dirac S(fij) function in the field strength. 

It is not surprising that topological A-B interactions are at the heart of quantum statis- 
tics. In its original form, the A-B effect [6] consists in the phase shift in electron interfer- 
ence due to the electromagnetic field, determined by the phase factor exp[(ze/ he) J A^dx 11 ] 
along a closed curve 7 passing through the beam along which the field strength vanishes. 
This effectj]] is counter-intuitive to the usual understanding that the influence of a classi- 
cal electromagnetic field on a charged particle can only occur through the local action of 
the field strength. In the context of quantum statistics, it means that non trivial statis- 
tics arise through topological "infinite" -distance interactions where no classical forces are 
present, as it should and as it is the case for Bose and Fermi statistics. Finally, singu- 
lar magnetic fields give an orientation to the plane, which, as already said, shows up in 

All this can be equivalently restated in a Lagrange formulation which describes again 
the system in topological terms, i.e. free particles minimally coupled to a potential vector 
whose dynamics is not Coulomb-like (Maxwell Lagrangian) but rather Chern-Simons [H] 



N 



N 1 r 

]T(-m^ + e(A(fi)^ - A°(fi))) - -e, up / d 2 rA^ A" (11) 
i=i " ^ 



with e^yp the completely antisymmetric tensor (the metric is (+,—,—), x M = (t, r) = 
(t,x,y), A^ = (A , A x , A y ), €ou = e 012 = +1). Solving the Euler-Lagrange equations, in 
particular 

^wkrm-^^p^ (12) 

1 The effect was first experimentally confirmed by R. G. Chambers [7], then by A. Tonomura [Hj. 
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leads to a magnetic field proportional to the density of particles in accordance with (10). 
Solving this last equation for A(f) in the Coulomb gauge gives 

Mr) = — E ^feff (13) 

j=l V 3 1 

in accordance with the A-B potential vector Q. Here again there is no Lorentz force, the 
potential vector is a pure gauge, the Chern-Simons term is metric independent, and the 
field strength is directly related to the matter current. 

Coming back to the Hamiltonian formulation ([8]), one might ask how the exclusion 
of the diagonal of the configuration space materializes in the Hamiltonian formulation. 
One way to look at it is perturbation theory [TOj [11]. Let us simplify the problem by 
considering the standard A-B problem, i.e. a charged particle in the plane coupled to a 
flux tube at the origin with the Hamiltonian 

Let us see what happens close to Bose statistics when a ~ (by periodicity a can always 
been chosen in [—1/2, +1/2], an interval of length 1, since in the one-body case one 
quantum of flux can always be gauged away via a regular gauge transformation). The 
A-B spectrum [6] is given by the Bessel functions 

iP(r) = e lW J\i- a \(kr) E = (15) 

with wavefunctions vanishing close to the origin r — > as J\i-. a \(kr) ~ r'' - "'. When the 
angular momentum 1^0, this is the only possible locally square-integrable function. 
However, when / = 0, one could have as well J_i a |(Ax) as a solution, since it is still 
locally square-integrable even though it diverges at the origin as r~' a L In principle, the 
general solution in the / = sector should be a linear combination of J\ a \ (kr) and J-\ a \ (kr), 
introducing an additional scale in the coefficient of the linear combination [12J. Restricting 



the space of solutions as in (15), i.e. wavefunctions vanishing at the origin, means that a 
short-range repulsive prescription has been imposed on the behavior of the wavefunctions 
when the particle comes close to the flux tube. One can give a more precise formulation 



of this fact by trying to compute in perturbation theory the spectrum (15). Expanding 



the square in the Hamiltonian (14), one finds that the a 2 /r 2 term, which is as singular as 
the kinetic term, is divergent at second order in perturbation theory in the / = sector. 
It follows that perturbation theory is not well defined in the problem as defined by the 



Hamiltonian (14). A renormalization has to be implemented: one realizes that by adding 



the countertem 7r|a|<5(f) to (14), i.e. by considering 



1 t-> kAf. 2 2-K\a\ . 
H=—ip-a— r ) 2 + —^-S{r) 16 
2m r z m 

the perturbative divergences due to the a 2 /r 2 term are exactly cancelled by those arising 
from the 7r|a|5(f) term at all orders in perturbation theory, giving back the spectrum 
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(15). Physically, this repulsive 5 contact term means that the particle is prevented from 
penetrating the core of the flux tube where the field strength is infinite, thus the (at 
least) r' a ' behavior when r — > 0. Note that this has been achieved without introducing 
any additional scale in the problem. 

Clearly, in the A^-body A-B anyon formulation of the model, the corresponding renor- 
malized Hamiltonian should read 

N 



1=1 j^i U i^j 



realizing the quantum mechanical exclusion of the diagonal of the configuration space in 
terms of contact repulsive interaction between particles. Note that the term 7r|a| T^i^j ^(^ij) 



in (17) can also be viewed [llj as the Pauli spin coupling of the spin of the particles to 



the singular magnetic field (10) associated to the flux tubes 



The anyon model defined in (17) is properly defined as far as short-distance consid- 
erations are concerned. It is the interacting formulation for regular wavefunctions of the 
free particles formulation for multivalued wavefunctions. Both Hamiltonians and H' N 
are equivalent, the former being more familiar in terms of usual quantum mechanics, the 
latter more relevant to study braiding and winding properties. 

The anyon model has been the subject of numerous studies in the eighties and the 
nineties [13] , some of them analytical, starting with the 2-body case which is solvable since 



its relative 2-body problem is the usual A-B problem (14) with an even (Bose) angular 
momentum I. The exact solution [TJ |5] for the relative 2-body problem is given by (15), 
I being an even integer, therefore when a is odd, I — a is odd, corresponding to Fermi 
statistics (the periodicity a — ► a±2p is manifest in the shift I — > I — a). These studies were 
followed by the 3-body [14] and then the A^-body problem [IS]. Statistical mechanics was 
also considered (second virial coefficient [T6l ITT], third virial coefficient jT8J). However, 
it soon became apparent that a complete A^-body spectrum was out of reach, to the 
exception of particular classes of exact eigenstates generalizing the 2-body eigenstates. 
Numerical [19] as well as semi-classical [20] studies were performed giving indications on 
the low energy A^-body spectrum. A systematic study of the model was achieved at first 
[2"T] and at second order [22J in perturbation theory (at second order the complexity of the 
model shows up clearly). Numerical studies [23] . taking some input from the perturbative 
results, were performed for the 3rd and 4th virial coefficients. Last but not least, on 
the experimental side, Laughlin quasiparticles [24] were put forward as the elementary 
excitations of highly-correlated fractional quantum Hall electron fluids [25J. They were 
supposed to carry a fractional charge and to obey anyon statistics [2H], a fact confirmed 
by Berry phase calculations, at least for quasiholes [27] (for quasiparticles the situation is 
less clear). The quasiparticles can propagate quantum-coherently in chiral edge channels, 
and constructively or destructively interfere. Unlike electrons, the interference condition 
for Laughlin quasiparticles has a non-vanishing statistical contribution which might be 
observed experimentally [25] . 

Some kind of simplification had to be made to render the model more tractable, and 
possibly solvable, at least in a certain sector. One realized that this was the case if one 
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considered, in addition to the singular statistical magnetic field, an external homogeneous 
magnetic field perpendicular to the plane, to which the charge of the anyons couple. In 
the case of a strong magnetic field, by projecting the system of anyons coupled to the 
magnetic field in its LLL, the model becomes solvable meaning that one can find a class of 
A^-body eigenstates which interpolates continuously from the LLL-Bose to the LLL-Fermi 
eigenstates basis: this is the LLL-anyon model [29] . 



2 The LLL-anyon model 

From now on, let us set the mass of the particles m = 1 and choose the statistical 
parameter a G [—2,0]. It is understood that all the results below are obtained for a 
in this interval, but they can be periodically continued to the whole real axis. Before 



introducing an external magnetic field, let us come back to the anyon Hamiltonian (17) 
and take advantage of wavefunctions vanishing at least as r~- a when — > (exclusion 
of the diagonal of the configuration space in the quantum mechanical formulation) by 
encoding this short distance behavior in the iV-body bosonic wave function [10] 



^(ri,r 2 ,...,r7v) = JJ^/^^i, f 2 , f N ) 



i<j 



iMri, f 2 , ■■■,tn) is regular but does not have to vanish at coinciding points. From in 
(17) one can compute the new Hamiltonian H N acting on ip(fx, f 2 , f N ). Since H N is 
itself obtained from the free Hamiltonian H' N in ^ via the singular gauge transformation 
([6]), it is more transparent to start directly from the free formulation. In comp 



(the free Hamiltonian is H' N — 
combined take the simple form 



^ Ei=i ^i) ^e wavefunction redefinitions (61) and (18) 



ex notation 



ip'{zi, z 2 , z N ; zi, z 2 , z N ) = Y\ Zij a 4>(zi, Z2, z N ; zi, z 2 , z N ) 



(19) 



i<j 



i<j Z ij 



in (19) encodes in the wavefunction the essence of 



The Jastrow-like prefactor f| 
anyon statistics: topological braiding phase and short-distance contact exclusion behavior. 
It is immediate that Hn rewrites as 



N 1 
H N = -2J2 did* + 2aJ2 $ - Bj) 

■ 1 ■ , ■ Zj 



(20) 



It is a non-Hermitian Hamiltonian (the transformation (19) is non-unitary), but it has 
a simple form, linear in a and well defined in perturbation theory (it is perturbatively 
divergence free). Any analytic wavefunction of the Zi S IS cl A^-body eigenstate of H N , 
and therefore of the A^-anyon Hamiltonian (17) taking into account (18). Analytical 



eigenstates are known to live in the LLL of a magnetic field, if such a field were present. Let 
us couple the electric charge of each anyon to an external magnetic field B perpendicular 
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to the plane such that by convention eB > and let us denote by u c = eB/2 half its 
cyclotron frequency. One now starts from the Landau Hamiltonian 



(21) 



In a magnetic field, the 1-body eigenstates have a long-distance Landau exponential be- 
havior exp(— l;UJ c ZiZi). Let us also encode this behavior in the wavefuction redefinition 
(19) so that it becomes 



N 



z 2 , z N ; zi,z 2 , z N ) = JJ exp(-^w c ^ z { z^{z x , z 2 , z N ; z l7 z 2 , 

i<j 



ZN 



i=l 



One obtains 



N 1 
H N = -2 Y^idA - uj c zA) + 2aJ2 ~ ^) + ^ 



1 = 1 



(22) 



(23) 



where the trivial constant energy shift from the Pauli coupling to the magnetic field has 
been ignored. As announced, acts trivially on iV-body eigenstates made of sym- 
metrized products of analytic 1-body LLL eigenstates 



(24) 



(in (24) the Landau exponential term is missing since it has already been taken into 
account in (22)). So, up to an overall normalization, 

N 



Tp(z 1 ,z 2 ,...,z N ;z 1 ,z 2 ,...,z N ) = Sym JJzj'; < h < l 2 < ... < l N (25) 

i=i 

is an eigenstate with a degenerate iV-body energy, En = Nuj c , a mere reflection of the 



fact that there are N particles in the LLL. From (22) and (25) one finally gets 

N N 



^ N N 

ip'izi, z 2 , z N ; z u z 2 , z N ) = Y[z~ j a exp(--uj c ^2z i z i )Sjm JJ^ 8 ; < h < l 2 < ... < l N 



8=1 



1=1 



(26) 

The basis (26) continuously interpolates when a = — > — 1 from the complete LLL-Bose 
iV-body basis to the complete LLL- Fermi iV-body basis. Indeed, when a = — 1, 

- N N 

i>\z x , z 2 , z N ; zi, z 2 , z N ) = exp(--u; c %Sym z l }\ < h < l 2 < ... < I 

(27) 



A' 



i=l 



i<j 



i=l 
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is equivalent to 



1 N N 
ip'(z 1 ,z 2 ,...,z N ;z 1 ,z 2 ,...,z N ) = exp(--u c ^2 z^) Antisym JJ^'; < l[ < 1' 2 < ... < l' N 

i=l i=l 

(28) 

i.e. the LLL fermionic basis. One has therefore obtained a complete LLL-Bose — >• LLL- 
Fermi interpolating basis which allows, in principle, for a complete knowledge of the 
LLL-anyon system with statistics intermediate between Bose and Fermi statistics. 

One could ask about going beyond the Fermi point a = —1 up to the Bose point 
a = —2. This question is related to the validity of the LLL projection, since ignoring 
higher Landau levels amounts to assuming that excited non LLL states above the iV-body 
LLL ground state have a non vanishing gap. Considerations around the Fermi point, as 
well as numerical and semiclassical analysis, support [29] this scheme as long as a does not 
come close to —2. However, when a — > —2, known linear as well as unknown nonlinear 
non LLL eigenstates do join the LLL ground state [31] • Said differently, the LLL-anyon 
basis (26) does not constitute a complete LLL-Bose basis when a — > —2, i.e. some N- 
body LLL bosonic quantum numbers are missing at this point. We will come back to this 
issue later. 



NlOr 




Figure 3: Linear and non linear non LLL eigenstates merge in the LLL ground state at 
the bosonic values of a. 



One has not seen yet any a dependence in the iV-body energy, a situation already 
encountered in the 1-body A-B problem, where the free continuous energy spectrum 
(15) is a- independent. This is due to the fact that a magnetic field does not confine 



particles: classical orbits are circular cyclotron orbits, but their centers, due to translation 
invariance, are located anywhere in the plane. Translation invariance in turn gives, in 
quantum mechanics, a Landau spectrum which is ^ independent, and therefore infinitely 
degenerate^] The degeneracy factor scales as the infinite surface V of the 2d sample: it 



Prom this point of view one can argue that the Landau spectrum is continuous, albeit being made of 
discrete Landau levels, due to the infinite degeneracy on each level. 
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is the flux of the magnetic field counted in units of the flux quantum <p = 2n/e (in units 
h=l) 

Nl = ^- (29) 

Statistical interactions being topological interactions, one does not expect, in the infinite 
plane limit, any effect on the iV-body energies. To see such an effect, one has to introduce 
a long-distance confinement, like putting the particles in a box. Let us rather introduce 
[30] a more convenient harmonic well confinement where the particles are trapped, so that 



the Landau Hamiltonian (21) becomes 

JV 



N R R N 

h' n = -2 j> - + + ^ J2 



(30) 



i=i 



i=l 



The virtue of the harmonic confinement is to lift the degeneracy with respect to the 
angular momentum ^ of the 1-body Landau eigenstates: the harmonic LLL spectrum^ 
becomes 



'a;, 



k+i 



71 



h>0; E = (uj t - uj c ){h + 1) +uj c 



(32) 



where uj t = \/uj 2 + tu 2 . Each harmonic LLL level in (32) has now a finite degeneracy, with 
an eigenstate still analytic in z^, up to the long-distance harmonic Landau combined ex- 
ponential behavior. Let us take into account this exponential behavior in the redefinition 
of the free iV-body wavefunction so that ( 22 ) now becomes 

1 N 

ip'(zi, z 2 , z N ; zi, z 2 , . ., z N ) = z tj a exp(--uj t ^ ZiZi)${zi, z 2 , ■■■ 



i<j 



i=l 



,z N ;z 1 ,z 2 ,...,z N ) 
(33) 



Starting from the Hamiltonian (30) one obtains [32| 133] 



N 



U t + LJ c _ - _ UJt - u c 



d i ) + 2ay2[-^—(d i -d j )- U ^-^]+Nu c 
^ l Zj — z- 9 J 



i<j 



(34) 

Again let us act on iV-body eigenstates made, in analogy with (25), of symmetrized 
products of the 1-body harmonic LLL eigenstates (32) 

TV 

V>(*i, z 2 ,...,z N ; zi_,z 2 , ...,z N ) = Sym JJ^'; < h < l 2 < ... < l N (35) 

i=i 

3 The complete 2d harmonic Landau spectrum is, with the convention eB > 0, 

Lo t (2m + k + 1) - liLu c ; ni>0, U e Z (31) 



The LLL quantum numbers are rij = and Zj > 
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Acting on this basis, the Hamiltonian (34) rewrites as 



N 

H N = (u t - u c ) [ J2 z i d i - <* 9 ~ + N ] + Nu c (36) 



so that the iV-anyon energy spectrum is 

N 



E N = (u t - u c ) [ k - u N{ - N l) +N}+ Nu c (37) 



i=i 



The iV-anyon spectrum (37) is a sum of 1-body harmonic LLL spectra shifted by the 
2-body statistical term — (uo t — u c )aN(N — l)/2. The effect of the harmonic well has 
been not only to lift the degeneracy with respect to the Zj's, but also to make the energy 
dependence on a explicit. When computing thermodynamical quantities like the equation 
of state, the harmonic well regulator will also be needed to compute finite quantities in 
a finite "harmonic" box, and then take the thermodynamic limit, by letting uj — > in an 
appropriate way 



The resulting eigenstates from (33) 



N N 



4>'(z 1 ,z 2 ,...,z N ;z 1 ,z 2 ,...,z N ) = J] Zjj 01 exp(-^ ^ z^)Sym J\z^; < h < h < ■■■ < In 



i<j i=l 



(ss; 

are called "linear states" since their energy (37) varies linearly with a. As already stressed 



they constitute a set of exact iV-body eigenstates which is only a small part of the complete 
iV-body spectrum, which remains mostly unknown. However, what makes, in the LLL 
context, these linear states particularly interesting is that they continuously interpolate 
when a = —>■ — 1 from the complete harmonic LLL-Bose basis to the complete harmonic 
LLL-Fermi basis. 

Before turning to LLL-anyon thermodynamics, let us reconsider the physical charge- 
flux composite interpretation of the anyon model, where the charges are now coupled to an 
external magnetic field. A given particle, say the iVth, sees a "positive" (eB > 0) magnetic 
field perpendicular to the plane, and N — 1 "negative" (e0 = 2na < 0, a G [—1, 0]) point 
vortices piercing the plane at the positions of the other particles. This is a screening 
regime: in the large iV limit where a mean field picture is expected to be valid, the more 
a is close to the fermionic point a = — 1, the more the external magnetic field is screened 
by the mean magnetic field associated with the vortices. In terms of the total (external + 
mean) magnetic field (B) that the iVth particle sees, or rather in terms of its flux V(B), 
or, when counted in units of the flux quantum, in terms of the Landau degeneracy (Nl), 
one has 

V(B)/ ( j )0 = (VB)/cj )0 + (N- l)0/0o i.e. (N L ) = N L + (N - l)a (39) 

Moving away from the Bose point, i.e. a < 0, as N increases the number (Nl) of 1-body 
quantum states available for the iVth particle in the LLL of (B) decreases. This sounds 
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reasonable, bearing in mind that a fermion occupies a quantum state to the exclusion of 
others (Pauli exclusion), whereas bosons can condense (Bose condensation). Introducing 
the LLL filling factor 

" = £ (40) 



one deduces from (39) a maximal critical filling [29] for which the screening is total, 
(N L ) = 

" = -- (41) 

a 

This is nothing but recognizing once more that bosons (a = 0) can infinitely fill a quantum 
state [y = oo), whereas fermions (a = — 1) are at most one per quantum state {y = 1). 
In between, one finds that there are at most — I /a anyons per quantum state. 

Interestingly enough, Haldane/exclusion statistics definition^] happens to coincide with 



(39): for a gas of particles obeying Haldane/exclusion statistics [M] with statistical param- 
eter g e [0, 1], given Nl degenerate energy levels and N — 1 particles already populating 
the levels, the number of quantum states still available for the iVth particle is given 



by (39) where —a is replaced by g 

d N = N L -(N- l)g (42) 



On the one hand, Haldane's definition (42) stems from an arbitrary combinatorial point 



of view, inspired by the Bose and Fermi counting of states. On the other hand, in the 



LLL-anyon model, (39) is obtained from a somehow ad- hoc mean field ansatz. We will 



come back to these issues in the next section. 

3 LLL-anyon thermodynamics 

Let us rewrite the iV-body energy ( J37| ) as [37] 

N 

E N = J2(e + kuj)-a ( % > ®<h<h<-<l N (43) 
i=i 

with uj = {uj t — u> c ) and e = oj c . Introducing the fugacity z and the inverse temperature 
(3, one wants to compute the thermodynamic potential 

oo 

lnZ(p,z)=\n(J2z N ZN)] Z Q = 1 (44) 

N=0 

where Z(/3, z) is the grand partition function defined in terms of the iV-body partition 
functions Zn = Tr exp(— (3H' N ) = Tr exp(— /3Hjf) = Tr exp(— (3H^). The thermodynamic 
potential rewrites as hiZ(fl,z) = Y^=i^>nZ n where, at order z n , the cluster coefficient 
b n only requires the knowledge of the Z^s, with % < n. One is interested in evaluating 

4 This is Haldane's statistics for one particle species. It can be generalized to the multispecies case. 
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the thermodynamic potential in the thermodynamical limit, i.e. to is small, which means, 
here, that the dimensionless quantity f3uj is small. The iV-body spectrum, as given in 
(43), allows to compute, at leading order in (3u —>■ 0, the Z^s for i < n, and thus the 6 n 's 

1 



1 e~ n ^ Vr 1 k + na 



@uj n 2 



n 

k=l 



k 



(45) 



One has still to give a meaning, in the thermodynamic limit (3uj = 0, to the scaling factor 
1/((3Q) in (45). To this purpose, one temporarily switches off the anyonic interaction and 
the external magnetic field, and considers a quantum gas of non interacting harmonic 
oscillators per se. One asks, when (3uj —>■ 0, for its cluster coefficients to yield the infinite 
box (plane wave) cluster coefficients. At order n in the cluster expansion, in d dimensions, 
one obtains [10J 

1 \l - V 



lim 



where A = ^2tt(3 is the thermal wavelength and V is the <i-dimensional infinite volume 
(in d — 2 dimensions, V is, as defined above, the infinite area of the 2d sample). Using 
the thermodynamic limit prescription (46), the cluster coefficient (45) rewrites, in the 
thermodynamic limit, as [29] 



b n = N L 



e -n(3u c «_1 k + na 



n 



n 

k=l 



k 



h = N L e~^ 



(47) 



The cluster expansion \nZ(P,z) = Y^=i^nZ n 1 as a power series of ze l3uJc < 1, can be 
summed up 

lnZ((3,z) = N L lny{ze-^ c ) 
where y(ze~ l3uJc ), a function of the variable ze - ^, is such that 



(46 



\ny = ze~^ c + ^ 



[ze 



n=2 



n 



n 

k=l 



It obeys 

and has in turn a power series expansion [38 



y - ze-^y 1 ^ = 1 



y 



ze 



n=2 



ze 



k + na 



k=2 



k 



(49) 



(50) 



(51) 



From (48) one infers that Z(/3,z) = y L so that 

oo 

Z((3, z )=y N ^ = l + N L ze~^ + N L ^ 



N 



ze 



-fko^N J-J k + N l + Na - 1 (52) 



N=2 



k=2 



k 
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Clearly, from (52), the iV-body partition function Zn is 



-N/3ui c 



N 

n 

k=2 



k + N L + Na - 1 
k 



(53) 



It is, by construction, positive. Necessarily, a and Nl being given, N has to be such 
that Nl + Na > 0. This always is the long as N is finite, since Nl scales like 

the infinite surface of the 2d sample. In the thermodynamic limit, where iV — > oo, the 
condition Nl + Na > implies for the filling factor 



v < 



1 



a 



(54) 



It is rather striking that the RHS of (54), which has just been derived from the exact 



computation of the cluster coefficients from the iV-body spectrum, is nothing but the 



critical filling (41) obtained in the mean field approach when the screening is total. 

The "degeneracy" associated with N anyons populating the LLL quantum states is, 



from (53) 



N, 



N 

n 

k=2 



k + Nl + Na — 1 N L (N + N L + Na - 1) 



k 



N\ 



(N L + Na)\ 



(55) 



-p + 



where a factorial with a negative argument has to be understood as (— p)\ = lim x _ 
x)\. 

When a = 0, this is the usual Bose counting factor for the number of ways to put N 
bosons in Nl states 

(N + N L -1)\ 

N\{N L -1)\ [ } 

When a = — 1, this is the Fermi counting factor Nl\/(N\(Nl — N)\). If one considers for 
a moment the statistical parameter to be a negative integer a < —1 , the degeneracy (55) 
still allows for a combinatorial interpretation [38J : provided again that Nl + Na > 0, 
it is the number of ways to put N particles on a circle consisting of Nl quantum states 
such that there are at least —a — 1 empty states in between two occupied states. When 
a = — 1, this is nothing but the usual exclusion mechanism for fermions (one fermion at 
most per quantum state). When a < — 1, i.e. beyond the Fermi point, more and more 
states are excluded between two filled states. In the case of interest a in [—1,0], one 
has a "fractional" exclusion where one can put more than one particle per quantum state 
according to the fractional a, but not infinitely many as in the Bose case. 

The degeneracy (55) originates from the exact iV-body spectrum (37). In the case of 



Haldane statistics as defined in (42), there is no Hamiltonian and no iV-body spectrum 



to begin with. One rather starts from the Bose counting factor (56) and bluntly replaces 



in accordance with (42), N L by N L — (N — l)g to obtain 



(N L -(N-l)(g-l))\ 
M(N L -(N- 1)0-1)! ' 



(57) 
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which indeed interpolates, when g = 1, to the Fermi counting factor. The degeneracy (57) 



is similar to (55): if one allows the exclusion parameter g to be an integer, it counts [38 



the number of ways to put iV particles on a line of finite length consisting of Nl quantum 
states such that there are at least g — 1 empty states in between two occupied states. Up 
to boundary conditions on the space of available quantum states (periodic versus infinite 



wall), both counting (55 57) are identical. In the thermodynamic limit when N becomes 
large, boundary conditions should not play a role anymore: not surprisingly, starting from 
(57) and following the usual route of statistical mechanics [ID] (saddle-point approxima- 



tion) leads, in the thermodynamic limit, to the same LLL-anyon thermodynamic potential 



given by the equations (48) and (50), where the anyonic parameter —a is replaced by the 
exclusion parameter g. 

Note that the grand partition factorization Z(j3,z) = y Nh in (48) could suggest 



an interpretation of y as a LLL-anyon grand-partition function for a single quantum state 
at energy u c , on the same footing as, when a = or a = —1, y = (1 =F ze~ l3uJc ) Tl is 
indeed the single quantum state grand partition function for a Bose or Fermi gas. This 
interpretation is not possible for the reason advocated above: it would yield, as soon as a 
is fractional, negative iV-body partition functions. This is clearly impossible: the iV-body 
anyonic system is, except in the Bose and Fermi cases, truly interacting and therefore its 
statistical mechanics is by no means factorisable to a single-state statistical mechanics. 

From (48, 50), the average energy E = —din Z(f3, z)/dj3 and the average particle num- 
ber N = zdln Z(/3, z)/dz or, equivalently, the filling factor v = N/Nl, can be computed. 
v satisfies 



y = l + 



v 



1 + av 



or, equivalently, using (50) 



ze 



(l + (l + a)v) 1+a (l + ais) 



(5? 



(59) 



When cr^O and a ^ 
in special cases like a - 



•1, this equation cannot in general be solved analytically, except 
— 1/2 (semions). The equation of state follows 



/3PV = ln(l 



1 + av' 



(60) 



In all these equations, it is understood from (54) that v < —1/a. When v = — 1/a, the 



pressure diverges, a manifestation of the fact that there are as many anyons as possible 
in the LLL, higher Landau levels being forbidden by construction. One also notes that, 



for the degenerate LLL gas, the filling factor in ( 59 ) is nothing but the mean occupation 
number n at energy e = uo c and fugacity z. As expected, (59) at a 



Bose mean occupation number n = ze 
Fermi mean occupation number n = ze 



-f3e 



/(I 



■ ze 
ze 



gives the standard 
whereas at a = —1 it gives the 



The entropy S = In Z((3, z) + (3E- (In z)N is (trivially E 
are in the LLL) 



Nu c since the N particles 



S = N L [(1 + v{\ + a)) ln(l + v(l + a)) - (1 + va) ln(l + va) - v\n.v\ (61) 
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It vanishes when v = — 1/a, an indication that the iV-body LLL anyon eigenstate is not 
degenerate at the critical filling. From (37), one infers that the iV-body eigenstate of 
lowest energy has all its one-body orbital momenta quantum numbers Zj = 0. It follows 



from (26) that, in the thermodynamic limit at the critical filling, the LLL-anyon non- 



degenerate groundstate wavefunction is 

1 N 

ip'(z 1 ,z 2 ,...,z N ;z 1 ,z 2 ,...,z N ) = Y[z ij a exp(--uj c ^2z i z i ); v 
with total angular momentum 



i<j 



i=l 



1 

a 



N(N — 1) 
2~v 



(62) 



(63) 



The pattern in (62) is reminiscent of the Laughlin wavefunctions at FQHE fillings v 
l/(2m+l) 



ip(zt, z 2 , z N ; zi, z 2 , ...,z N ) = J j 



4 m+lex p( 



1 N 



1 



i<j 



i=l 



2m + 1 



(64) 



On the one hand, Laughlin wavefunctions are fermionic, their filling factors are rational 
numbers smaller than 1, and they are approximate solutions to the underlying iV-body 
Coulomb dynamics in a strong magnetic field. On the other hand, LLL-anyon wavefunc- 
tions are multivalued, their filling factor continuously interpolates between oo and 1, and 
they are exact solutions to the iV-body LLL anyon problem. Still, the similarity between 



(62) and (64) is striking. 



Trying to push (62 ) further beyond the Fermi point eventually up to the Bose point at 
a = —2, one obtains a Bose gas at filling v = 1/2 with the non-degenerate wavefunction 



V>'(zi, z 2 , z N ; Zx,z 2 , z N ) = Yl 



4-exp( 



1 N 



V 



(65) 



i<j 



i=l 



One already knows that the LLL-anyon basis (26) is not interpolating to the complete 
LLL-Bose basis when a = —2. At this point, non LLL iV-body eigenstates merge in the 
LLL ground state to compensate for some missing bosonic quantum numbers -see Figure [3j 
Clearly, (65) should reproduce, by periodicity, the bosonic non-degenerate wavefunction 



(62 ) at a = 0, but it does not. On the same footing, when a = —2 the critical filling should 
be bosonic, i.e. v = oo, whereas v — 1/2. The unphysical critical filling discontinuity, 
oo versus 1/2, is yet another manifestation of the missing bosonic quantum numbers. 
In other words, the very eigenstates which join the LLL ground state at the Bose point 
a = —2 and provide for the missing quantum numbers, have the effect to smooth out the 
critical filling discontinuity. Still, it has been shown [35] that the stronger the magnetic 
field B is, the more valid (62) remains closer and closer to a = —2. The limit a — > —2 
is, due to periodicity, the same as the limit a — > from above, which can be described as 
an anti-screening regime. One concludes that close to the Bose point a = 0, the critical 
filling of a LLL-anyon gas isz/ = ooorz/=l/2 depending on infinitesimally moving 
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away from the Bose point in the screening regime (the ground state wavefunction is the 
usual non degenerate bosonic wavefunction), or in the anti-screening regime (the ground 



state wavefunction is (65)). Again, the Bose point has a somehow singular behavior, a 



feature already encountered in perturbation theory Note finally that the occurrence of 
the v = 1/2 fraction for the bosonic filling factor in the antiscreening regime is physically 
challenging: fast rotating Bose-Einstein condensates in the FQHE regime are expected 



[36] to reach a 1/2 filling described by the Laughlin-like wavefunction (65). 




Figure 4: The critical LLL-anyon filling curve as a function of a. The critical Bose filling 
v = | occurs at the Bose points in the anti-screening regime. The continuity of the 
critical curve at these points is restaured by the non LLL eigenstates joining the LLL 
ground state. 



So far one has been concerned with two-dimensional systems: in the thermodynamic 
limit, a single particle in the LLL, and, consequently, a gas of LLL-anyons, are two dimen- 
sional, as can be seen from the Nl — V scalings^] of the 1-body LLL partition function 
Zlll = Niexp(—Pu c ) and the LLL anyon thermodynamic potential (48). Denoting by 
Plll(c) = NL5(e — u c ) the 1-body LLL density of states, (48) can be rewritten as 



hi Z (13, z) 



p LLL (e)\ny(ze 0e )de 



(66) 



Convincingly, in (66) the one-body dynamics of individual particles is described by the 



one-body density of states, whereas the LLL anyon statistical collective behavior is en- 
coded in the y function which depends on the statistical parameter a. 

One might ask about other integrable iV-body systems which would lead to the same 
kind of statistics. It would be tempting to define a model obeying fractional/exclusion 
statistics if, its one-body density of states p(e) being given, its thermodynamic potential 
has the form 

/■oo 

\nZ(f3,z)= p(e)\ny(ze~^)de (67) 



with 



y - ze-? e y l+a 



(68) 



3 In the LLL there is only one quantum number U per particle, still the system is 2d. 



19 



so that 

oo n , 

y = l + ze-^ + Y J ^r\{- ± P (69) 

n=2 k=2 

The mean occupation number follows as n = zdlny/dz. It obeys 

Tl V — 1 

V = l + — ; or n=—?- ->0 (70) 

1 + an 1 — a(y — 1) 



or, equivalently, 



ze -P* = - (71) 

(1 + (1 + a)n) 1+a n + an)-" 1 ; 



One has the duality relation j?T] 



y y 

or, equivalently 



1 = - + -; where y - (^"^"V" 1 "" = 1 (72) 



— an h = 1 (73) 

a 



where h is related to y as n to y in (70). The duality relation ( 72|73 ) can be interpreted 
as a particle-hole symmetry relation. Setting t = ze~ f3uJc , one also has a simple expression 
pE2] for dn(t) /dt 

dft 

t— = nil + (1 + a)n)(l + an) (74) 
at 

All these equations have been understood as arising microscopically from the LLL 
anyon Hamiltonian with one-body density of states p(e) = Plll(c). It happens that it 
is possible to find another iV-body microscopic Hamiltonian which leads to the thermo- 



dynamics (67). Consider, in one dimension, the integrable iV-body Calogero model 



with inverse-square 2-body interactions 



N o2 N 

n N 



where Xi represents the position of the i-th particle on the infinite Id line. This model 
is known to describe particles with nontrivial statistics in one dimension interpolating 
from Bose (a = 0) to Fermi (a = —1) statistics. It means that the 1/x 2 Calogero 
interaction is purely statistical, without any classical effect on particle motions, up to a 
overall reshuffling of the particles |Sj. The Calogero model remains integrable when, as 



in (75), a confining Id harmonic well is added. This is the harmonic Calogero model, 
whereas the Calogero-Sutherland model [16] would have the particles confined on a circle. 
The effect of the harmonic well is, as in the LLL anyon case, to lift the thermodynamic 
limit degeneracy in such a way that the iV-body harmonic Calogero spectrum ends up 
depending on the Calogero coupling constant a 

E N = uj[J2h-a N ^ N ~^ +y]; < h < l 2 < ... < l N (76) 
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Here the Zj's correspond to the quantum numbers of the 1-d harmonic Hermite polynomials 
free 1-body eigenstates 



UJXi 



k>0; E = u(U + 



(77) 



It is remarkable that (76) happens to be again of the form (43) with Co = lo, €q = to/2. 



Following the same steps as in the LLL-anyon case, and using again (46) while taking the 
thermodynamic limit (3uj — > 0, the Calogero cluster coefficients rewrite as 



K 



L 1 



n\/n 



n-l 

n 

fc=i 



k + na 



L 
A 



(78) 



where the infinite length of the Id line has been denoted by L. The cluster expansion can 



still be resumed using (49) provided the unwanted 1/y/n term in (78) is properly taken 



care of. Introducing the Id plane wave momentum k 



A 



n 



1 

2vr 



dke 



-n/3* 



(79) 



and denoting the 1-body energy as e = k 2 /2, one finally obtains 



\nZ{P,z) 



p (e)lny(ze pe )de 



(80) 



where 



Po(eJ 



7T' 



This is not a surprise: in the 
= kf/2 kept fixed, the Hermite 



is the free 1-body density of states in one dimension, 
thermodynamic limit uo — ► 0, where l{ — > oo with liio 
polynomial becomes a plane wave of momentum ki. 

From (80), one conclude^] that, in the thermodynamic limit, the Calogero model has 

151 according to (67) and (68), interpolating, 

-1. 



indeed a LLL-anyon/exclusion like statistics 
as it should, from a free bosonic Id gas at a 



to a free fermionic Id gas at a ■ 
It follows that the 2d LLL-anyon and Id Calogero models, which seem a priori unre- 
lated, do obey the same type of statistics. This is not a coincidence. Looking at their 



limit of the former. This remains true in the thermo- 



harmonic iV-body spectrum (37) and (76), one realizes that, up to an irrelevant zero-point 
energy, the latter is the B 
dynamic limit u 
B 



0. So, not only (66) and (80) are of the same type, but also, when 



0, (|66|) has to become (80). It follows that, necessarily, the 1-body densities of states 
satisfy lim B ^ p LLL (e) 



PLLL(e) and p { 



i.e. 



.. eBV.. eB, 

lim o(e 

b^o 2vr v 2 ' 



L 



-kV2~€ 



(82) 



6 The same conclusion would be reached starting form the Calogero-Sutherland model and taking the 
corresponding thermodynamic limit, i.e. the radius of the confining circle going to infinity. 
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a relation which has to be understood as arising in the thermodynamic limit to — > 0. 



To arrive at (82), one could as well consider directly the 1-body harmonic LLL spec- 



trum (32) and harmonic Id spectrum (77) 



E = (uj t - u c )(li + 1) + u) c \ E 



u{k + 2 ) 



33) 



They are such that the latter is the vanishing B limit of the former, so it is the case for the 
corresponding 1-body partition functions. Taking therQthe thermodynamic limit (3uj — > 0, 
i.e. (46), implies the relation riniB_>o Zlll = Z , where Zlll is, as above, the LLL 
partition function and Zq is the free partition function in one dimension. Consequently 



for the densities of states (the inverse Laplace transforms) the relation (82) follows. This 



result has its roots in the different energy gaps of the spectra (83) at small 10: in the 
harmonic LLL case, the gap behaves like u 2 /(2lj c ), whereas, in the Id harmonic case, the 
gap is us. 



The relation (82) could also have been understood from the 1-body eigenstates them- 
selves. In the limit B — > 0, the LLL induced harmonic analytic eigenstates are, from 



(32) 



h+i 



z/e 



34) 



There is only one parameter ui left so that the states in (84) can be put in one-to-one 



correspondence with the Hermite polynomials (77) via the Bargmann transform 



Voo k+1 z l ; 



35) 



From (85) one can infer [47] that the iV-body harmonic anyon eigenstates (38) are a 



coherent state representation of the iV-body harmonic Calogero eigenstates. 

From all these considerations (thermodynamics, eigenstates,...) it follows that the 
vanishing magnetic field limitj^] of the LLL-anyon model is the Calogero model itself. It 
seems paradoxical to consider such a limit in the LLL which assumes a priori a strong 
magnetic field. Still, doing so, one has dimensionally reduced the 2d anyon model to the 
Id Calogero model. This dimensional reduction has a simple geometrical interpretation. 



The LLL induced harmonic states (84) are localized in the vicinity of circles of radius U/uj. 
In the thermodynamic limit, one has U —>■ oo with l^u = kf/2 kept fixed. It follows that 
the corresponding Id Hermite polynomials Hi t , which become in this limit plane waves 
of momentum hi, have a radius of localization diverging like kf/u 2 . The dimensional 



The order of limits is crucial here: first the limit B — > 0, then the thermodynamic limit uj —> 0. 
8 Since one has ended up by taking the limit B — > 0, one could have avoided right from the beginning 
to introduce a B field, and started directly from the harmonic TV-body anyon model. What has been 
done above by taking the limit B — > is nothing but to project the harmonic anyon model on the LLL 



induced harmonic subspace ( 84 ) (the B field and its LLL should still be invoked to justify the selection 



of the LLL quantum numbers in the 2d harmonic basis) and to recognize that the projected harmonic 
anyon model is the harmonic Calogero model. This relation remains true in the thermodynamic limit 
uj -> 0. 
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reduction which has taken place consists in going at infinity on the edge of the plane: in 
the thermodynamic limit, the Calogero model can be viewed as the edge projection of the 
anyon model. 

The LLL anyon thermodynamics, or, equivalently, the Haldane/exclusion thermody- 
namics, and the Calogero thermodynamics as well, have been the subject of an intense 
activity since the mid-nineties. Let us mention their relevance in more abstract contexts, 
such as conformal field theories |48j . On the experimental side, FQHE edge currents can 
be modelled by quasiparticles with fractional statistics, which in turn might affect their 
transport properties such as the current shot noise [1HJ E2] • 



4 Minimal Difference Partitions and Trees 

Up to now one has been concerned with quantum mechanical models defined by a micro- 
scopic quantum Hamiltonian. Both the LLL anyon and Calogero models have been shown 
to have a thermodynamics controlled by (67) and (68). Let us leave quantum mechanics 
and address a pure combinatorial problem, the minimal difference partition problem [50] . 
Consider the number p(E, N) of partitions of an integer E into iV integer parts where 
each part differs from the next by at least an integer p and the smallest part is > I. Usual 
integer partitions correspond to p = and I = 1, whereas restricted partitions, where the 
parts have to be different, correspond to p = 1 and 1 = 1. 



h 



> p 



Wi 



>p 



>p 



> 



p 



> I 



Figure 5: A minimal difference partition configuration, or Young diagram. The column 
heights are such that (ij — > p for i = 1,2, . . . , N — 1 and h > I. Their total height 
is E = J2iLi h- Wh is the width of the Young diagram at height h, i.e. the number of 
columns whose heights > h. 
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It is known that 



x lN+pN(N-l)/2 



X 



[1 — x)(l — x 2 )...(l — X 



(86) 



The p(E,N) generating function Z(x,z) = N p(E , N)x E ] z factorizes when p = or 
p = 1 

OO -. oo 



p = 0, Z(x, z) 



i=0 



x l+l z ) 



7) 



i=0 



In terms of bosons or fermions, (87) is the grand partition function for a bosonic or 
fermionic gas with fugacity z and, denoting x = e -13 , temperature T = 1/(3 where 



71: 



i=0 



i=0 



with rii = 0,1,2,... in the Bose case (p = 0) and 7%i — 0,1 in the Fermi case (p = 1). 
Equivalent ly 



N 



e = y,u 



(89) 



i=i 



with I <h <h < ■ ■ ■ <In (Bose) or / < l x < h < ■ ■ ■ < In (Fermi) 



When p is an integer > 2, (86) can be regarded as the iV-body partition function of 



an interacting bosonic gas with the iV-body spectrum 

N 



E = ^2k+pN(N -l)/2; l<h<l 2 < ...<In 



(90) 



i=l 



Clearly, (90) goes beyond the Fermi point p — 1 and describes some kind of "super- 
fermions". In contrast to the Bose and Fermi factorization such as (|87|) is not 

on 

(91) 



possible, due to the interacting nature of (90). One has instead the functional relation 

Z(x, z) = Z(x, xz) + x zZ(x, x p z) 



which embodies the combinatorial identity 



p{E,N)=p (E-p 



N 



(92) 



where p (E,N) stands for the usual partition counting. 

One could push (51] this analysis further to p real positive. When p e [0, 1] and 1 = 1, 
one would obtain a partition problem interpolating between the usual (bosonic) one and 
the restricted (fermionic) one. It is manifest that, if p is replaced by —a, the spectrum 
(90) coincides, under a rescaling and up to an irrelevant zero-point energy, with the iV- 
body quantum spectrum (76) of the harmonic Calogero model. In a partition problem, 
one is interested in the large E and iV asymptotic behavior of p(E, N), which corresponds 
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to the regime x — > 1, i.e. (3 — > 0. Consider the cluster expansion lnZ(x, z) = Yl^Li bnZ r ' 
In the limit /3 — > one obtains 



J e -nZ/3 



n-1 

fe=i 







(93) 



The limit /3 — > should not be confused with the thermodynamic limit in quantum 
systems. There is no thermodynamic limit prescription like (46). Still, using (49) (with 
uj c replaced by I) and taking care of the unwanted 1/n factor in (93), one obtains, provided 
that ze-& < 1, 



\ny(ze-P £ )de 



with 



y - ze- pe y l - p 



1 



(94) 



(95) 



This is again of the form (67) and (68), the statistical parameter —a being replaced 
by the minimal difference partition parameter p, and the 1-body density of states being 
the Heaviside function p(e) = 8(e — I). The minimal difference partition combinatorics 
is equivalently described, in the small (3 limit, by a gas of particles obeying exclusion 
statistics with a uniform density of stateg^J 

This correspondence happens to be useful technically: (94) and (95) are the build- 
ing blocks of the minimal difference partition asymptotics. The average integer E = 
—d In Z(/3, z)/d(3 = nede and the average number of integer parts N = zd In Z((3, z) / dz 
f, nde, are both given in terms of n = zdlny/dz, the mean occupation number at "part" 
e and fugacity z, which satisfies 



ze 



-f3e 
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One obtains 



(1 + (1 - p)ny-p(l - pn)P 
1 



with n < 



1 

V 



E-IN = - In Z(P,z 
P 



N 
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\ny(ze 



(96) 



(97) 



so that the entrop> S = \nZ((3, z) + f3E — (In. z)N rewrites as 



S = 2f3(E-lN 



2 



-PNy 



(98) 



with 



E-IN - -N 2 
2 



1 

T 2 



ln(l - u) 



du 



u 



(99) 



Inverting (99) gives (3 as a function of E and N so that the entropy S in (98) becomes a 
function of E and N only. Doing so, one has a definite information [51] on the asymptotic 
behavior of p(E, N) ~ e s( * E ' N ^ when E and iV are large, and also, of p(E) = Y1'n=i p{^i N) 



9 There is no microscopic quantum Hamiltonian leading to (94 1 and (95). 
10 The simple expression in (97 1 for N is possible because of tne constant density of states. 
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when E is large. One obtains a generalization of the Hardy- Rama junan asymptotics [52] 
to the minimal difference partition problem. One can also obtain [53] the average limit 
shape of the Young diagrams associated with the minimal difference partition problem, 
generalizing the usual partition limit shape [M]- The limit shape at a part of height h 
depends solely on the statistical function y evaluated at e = h and at z = 1 



PW h = \ny(e 



(100) 



where (3 scales as (3 2 E = J °° \ny(e~ e )de. 

So far p being a positive integer has insured that the iV-body spectrum in (90) is well 
defined. However, y in (95) is still meaningful when p is a negative integer. It is the 



(1 — p)-aiy tree generating function, so that the coefficient at order n of its expansion in 
powers of ze~@ e as given in (69) (with —a replaced by p) is the number of ways to build a 
(1 — p)-ary tree with n nodes. For example, at p = — 1, y generates the Catalan numbers 
associated with binary trees. 
Consider, toy model [55] 



the factorized (1 — p)-ary tree generating function 



Z(x,z) = l[y(zx l+l ) 



(101) 



=o 



where y satisfies (101) with e 
case). Its combinatorial interpretation is that p(E,N) deduced from (101) counts the 



I + i. (101) narrows down to (87) when p = (Bose 



number of usual partitions of an integer E into iV integer parts bigger or equal to I, with 
an additional degeneracy stemming from the (1 — p)-tree arborescence when, in a given 
partition, a part occurs n times. This enlarged degeneracy goes beyond the Bose point to 
define some kind of " superbosons" . 

One can analytically continue p to the whole negative real axis. In the large E and 
N limit, i.e. f3 smaller and smaller, one encounters a maximal temperature beyond which 
it is not possible to heat the system. Indeed, from (95) it follows that y(zx l+t ) in (101) 



obeys to y — ze 



'V 



1, which is well defined only if 



ze 



-pi 



< (i-p) p -\-py p < i 



When z — 1, it defines a dimensionless "Hagedorn temperature" 

T 1 



[1 — p) ln(l — p) + p ln(— p) 



(102) 



(103) 



just below which E and N become large so that the asymptotic of p(E, N) can be ad- 
dressed. 



5 Conclusion 



In two dimensions intermediate anyonic statistics interpolating from Bose to Fermi statis- 
tics are allowed. Their definition does not involve anything else than the usual concept 
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at the basis of quantum statistics, namely free particles endowed with particular bound- 
ary exchange conditions on their A-body wavefunctions. It happens that these boundary 
conditions have a much richer structure in two dimensions than in three and higher dimen- 
sions. This in turn can be understood in terms of the topology of paths in the iV-particle 
configuration space, where non trivial braiding occurs in two dimensions, and not in higher 
dimensions. A flux-charge composite picture emerges to encode the braiding statistics in 
physical terms, via topological Aharonov-Bohm interactions and singular magnetic fields. 

The anyon model as such is certainly fascinating as far as quantum mechanics is con- 
cerned, but it remains an abstract construction whose complexity is daunting. However, 
when projected onto the LLL of an external magnetic field, the model becomes tractable 
and, even more, solvable. The LLL set up is clearly adapted to the QHE and to the FQHE 
physics. Haldane/exclusion statistics, which can be obtained as a LLL-anyon mean-field 
picture in the screening regime, leads to LLL-anyon thermodynamics. 

It would certainly be rewarding if LLL anyons could be relevant experimentally, for 
example by uncovering some experimental hints at FQHE filling factors of the existence of 
quasiparticles with anyonic/exclusion statistics. Fractional charges have already been seen 
in shot noise FQHE experiments [56], but the nontrivial statistical nature of the charge 
carriers in FQHE edge currents has so far remained elusive in experiments which rely 
mainly on Aharonov-Bohm interferometry [2H]- Note also a recent proposal for the pos- 
sible experimental tracking of abelian and nonabelian anyonic statistics in Mach-Zehnder 
interferometers |57j . 

Finally, on the theoretical side, physical interactions, together with topological anyonic 
interactions, should also be taken into account in order to produce more realistic models. 
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